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ABSTRACT 


^The  life  prediction  and  reliability  assessment  of  the  M60  tank  torsion 
bars  were  obtained  from  applying  a  methodology  involving  several  disciplines. 
The  disciplines  included:  structural  mechanics,  fracture  mechanics, 
statistical  reliability,  mechanical  testing,  nondestructive  examination 
(NDE),  quality  assurance  (QA),  and  metallurgical  and  fractographic 
evaluations. 

The  methodology  was  applied  in  order  to  obtain  a  procedure  for  increas¬ 
ing  torsion  bar  life  and  enhance  the  reliability  of  the  M60  tank.  U.S.  Army 
reliability  performance  is  commonly  measured  in  mean  miles  between  failures 
(MMBF) .  In  order  to  introduce  a  more  meaningful  measure  of  component 
acceptability,  this  report  describes  methods  for  obtaining  minimum  life 
estimate  at  a  specified  probability  of  survival  using  the  Monte  Carlo  pro¬ 
cess.  This  includes  predicting  remaining  life  with  known  risk  and  an 
approach  to  perform  trade-off  costs  or  redesign  for  increased  life.  The 
result  is  a  more  meaningful  measure  of  component  acceptability.  The 
reliability  versus  bar  life  computation  indicated  a  negligible  amount  of 
life  after  a  bar  flaw  was  initiated  in  the  bar.  The  minimum  life  estimates 
prior  to  crack  initiation  were  in  good  agreement  with  frequency  of  actual  bai 
failures  recorded  at  Aberdeen  Proving  Grounds.  i  o  i ;  <.<  1  t  _  ;  |u  > 
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I.  INTRODUCTION 

The  M60  tank  (Figure  1),  which  was  developed  in  the  1960's,  is  still  the  prin¬ 
cipal  Army  heavy  tank.  It  is  currently  undergoing  a  major  "Rebuild"  program  to 
increase  its  effectiveness.  The  M60  program  manager  has  expressed  concern  regarding 
the  reliability  of  various  components  of  the  overhauled  tanks  as  a  result  of  failures 
experienced  during  proving  ground  tests.  In  a  report  to  Congress,  January  29,  1980, 
the  GAO  expressed  a  similar  concern  on  the  reliability  of  other  systems.  This  study 
is  an  effort  by  an  interdisciplinary  task  group  assembled  at  AMMRC  from  specialists 
in  material  science,  engineering,  structural  mechanics,  nondestructive  testing,  and 
statistical  analysis.  The  charter  of  the  task  group  was  to  develop  a  methodology 
for  life  prediction  of  structural  components  in  the  M60  tank.  The  torsion  bar 
(Figures  2  and  3)  was  chosen  because  of  the  availability  of  a  data  base  obtained 
from  a  comprehensive  research  at  AMMRC,  proving  grounds  and  overhaul  depots.  It 
should  be  noted  that  the  torsion  bar  performs  like  a  spring  by  providing  the  neces¬ 
sary  motion  for  the  wheel  to  follow  the  terrain.  The  M60  has  six  torsion  bars  on 
each  side.  The  front  and  rear  wheel  arms  have  shock  absorbers  and  stops,  while  the 
intermediate  wheel  arms  have  only  stops  to  limit  their  motion.  Torsion  bar  failure 
can  lead  to  reduced  mobility  and  track  throwing,  or  failure  to  other  overloaded  bars 
and  eventual  total  loss  of  mobility  of  the  tank. 

The  task  groups  made  visits  to  manufacturing  facilities,  the  overhaul  depot, 
field  test  sites,  AMSAA,  and  TACOM  in  order  to  obtain  data  and  background  information. 
Detailed  failure  analysis  metallurgical  examination,  material  testing,  NDE,  detailed 
stress  analysis,  and  a  comprehensive  probablistic  analysis  were  completed  by  the 
group.  The  paper  will  discuss  in  detail  the  above  studies,  including  major  findings 
of  a  sensitivity  analysis,  and  the  final  conclusions  with  regards  to  design  improve¬ 
ments,  material  requirements,  and  the  feasibility  of  predictive  methodology  for  tor¬ 
sion  bars.  The  methodology  used  here  circumvents  the  present  deterministic  approaches 
used  in  establishing  design  life.  It  simulates  variability  of  loading  and  material 
by  introducing  Monte  Carlo  methods. I 

II.  RELIABILITY  AND  LIFE  PREDICTION  METHODOLOGY 

The  approach  taken  to  developing  life  prediction  methodology  was  basically  to 
obtain  as  much  data  about  the  torsion  bar  as  the  program  would  allow,  analyze  the 
data,  perform  the  reliability  study,  and  make  recommendations. 

The  M60  task  group  visited  various  manufacturers  of  the  torsion  bars  and  obtained 
quality  assurance  procedures,  test  results,  manufacturing  procedures,  and  machining 
methods  (General  Dynamics  and  FMC).  The  developer,  the  Army  Tank-Automotive  Command 
(TACOM),  provided  the  task  group  with  the  requirements  and  specifications,  Depot 
Maintenance  Work  Requirement  (DMWR),^  described  the  user  experiences.  Aberdeen  Proving 
Ground  (APG)  test  facilities  and  the  Materials  Testing  Directorate  provided  spectrum 
load  data  from  a  specific  course.  They  also  carefully  monitored  bar  failures.  The 
Anniston  Army  overhaul  depot  described  inspection  procedures  information,  NDE,  and 
replacement  requirements.  A  major  part  of  the  material  data  generated,  involved 
mechanical  and  various  metallurgical  properties,  and  their  relation  to  the  speci-  i 
fication.  These  properties  were  then  incorporated  with  the  data  bank  obtained  from 
the  QA  testing  of  the  manufacturer  and  published  data  on  the  material. 

1.  NEAL,  D.  M.,  and  MASON,  D.  S.  Determination  of  Structural  Reliability  Using  a  Flaw  Simulation  Scheme.  Army  Materia ti  and 
Mechanic*  Reiearch  Center,  AMMRC  TR  81-53,  October  1981. 

2.  Suspension  System  Components  for  M60.  Depot  Maintenance  Work  Requirement,  TACOM,  DMWR9-2350-253-1-5. 
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The  methodology  then  required  accurate  analysis  of  the  mechanical  loading,  stresses, 
fatigue,  and  fracture  evaluation  of  the  torsion  bar. 


Finally,  all  the  data  is  incorporated  in  a  Monte  Carlo  method  in  order  to 
determine  the  reliability  of  the  torsion  bar.  This  method  requires  modeling  of  the 
spectrum  loads  and  the  material  fatigue  life  with  respect  to  crack  growth  or  stress/ 
cycles  to  failure  (S/N).  The  reliability  of  a  single  torsion  bar  (or  system)  can 
then  be  studied  under  various  assumptions  of  terrain,  speed,  improvement  in  designer 
material,  inspection  procedure,  replacement  period  or  any  combinations. 


III.  MATERIAL  CHARACTERIZATION  DATA 

Several  sources  were  used  in  obtaining  material  properties.  These  covered  QA 
data  obtained  from  the  manufacturer,  mechanical  and  metallurgical  tests  performed  at 
AMMRC  by  the  task  group,  published  material  data,  and  discussions  with  the  other 
torsion  bars  manufacturers. 

The  M60  torsion  bars  are  currently  being  manufactured  by  General  Dynamics  (GD, 
Scranton,  PA)  and  Machine  Products  Co. ,  according  to  MIL-5-45387A,  from  AISI  8660 
steel.  It  should  be  noted  that  the  bar  specifications  do  not  have  an  explicit 
requirement  for  fracture  toughness  or  Charpy  impact  energy  and  the  only  requirement 
is  placed  on  Rockwell  Hardness.  The  metallurgical  data  reported  here  was  obtained 
from  the  failed  torsion  bars  obtained  from  Aberdeen  Proving  Grounds  (APG)  and 
Anniston  Army  Depot. 

Fatigue  Data 

As  part  of  the  quality  assurance  (QA)  requirements,  several  torsions  bars  from 
each  heat  were  tested.  The  torsional  test  machine  is  similar  to  that  shown  in 
Figure  4.  The  QA  data  obtained  from  GD,  shown  in  Table  1,  was  that  of  the  number  of 
cycles  to  failure  of  the  torsion  bar  for  a  specified  prescribed  angle  of  twist  (42°) 
which  is  based  on  the  maximum  allowable  torsional  shear  stress  of  140  ksi.  Table  1 
also  shows  the  cause  and  location  of  failure  in  each  test.  This  data  is  represented 
by  the  probability  of  failure  versus  number  of  cycles  to  failure  (Figures  5  through 
8)  for  the  individual  years  1978-1982.  The  Kruskal-Wallis  multi-sample  test  for 
identical  populations,-*  showed  that  there  are  no  significant  differences  among 
individual  yearly  fatigue  test  results,  therefore,  the  data  was  pooled  together. 

The  data  shown  in  Figure  9  is  the  accumulation  from  test  results  on  the  various 
heats  of  production  from  1977-1982.  The  1983  data  were  not  included  as  they  represent 
a  new  bar  modification. 

The  QA  specification  requires  that  three  bars  be  tested  for  every  heat  with  at 
least  45,000  fatigue  cycles  obtained  prior  to  bar  failure.  The  approach  currently 
used  by  the  manufacturers  is  to  run  the  test  until  the  torsion  bar  fails  (Table  1). 

In  many  instances  the  fatigue  test  machine  has  failed,  due  to  the  excessively  large 
number  cycles  required  resulting  in  an  increased  cost  for  quality  control.  In  order 
to  overcome  this  difficulty  a  statistical  procedure  was  applied  using  a  Weibull 
censored  data  analysis^  approach.  This  procedure  obtains  the  design  allowable 


3.  HOLLANDER,  M.,  and  WOLFE,  D.  A.  Nonparametric  Statistical  Methods.  John  Wiley  A  Sons,  New  York,  1973. 

4.  LAWLESS,  J.  F.  Contruction  of  Tolerance  Bounds  for  the  Extreme  Value  and  Weibull  Distributions.  Technometrics,  v.  17,  no.  2, 
May  1975. 
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values  for  failure  data  plus  run-out  values.  The  run-out  values  are  non-failed  bars 
with  only  the  information  that  they  have  not  failed  at  a  specified  number  of  cycles. 

Using  the  QA  data  from  GD  it  was  determined  that  the  maximum  number  of  required 
cycles  can  be  reduced  to  110,000  as  compared  to  the  maximum  of  570,000  currently 
being  required  to  cause  failure  of  the  most  durable  bar.  A  comparison  of  material 
design  allowable  values  obtained  from  both  complete  sample  (all  data)  and  censored 
data  (only  failed  data  below  110,000  cycles)  showed  excellent  agreement.  The  allow¬ 
able  represents  a  value  determined  from  a  specified  probability  of  survival  with  a 
95%  confidence  in  the  assertion.  Survival  probabilities  are  0.99  for  the  "A"  allow¬ 
able  and  0.90  for  the  "B"  allowable. 

The  above  indicates  that  equivalent  information  could  be  obtained  with  a  large 
saving  in  time  while  reducing  fatigue  machine  failures.  Figures  10  and  11  show  the 
censored  data  results. 

Metallurgical  Evaluation  and  Other  Mechanical  Properties 

A  metallurgical  evaluation  was  conducted  on  two  failed  M60  torsion  bars,  one 
produced  by  General  Dynamics  (GD)  and  the  other  by  Machine  Products  Company  (MPU). 
The  parameters  addressed  in  this  evaluation  are: 

o  chemical  analysis 
o  light  microscopy 
o  retained  austenite 
o  hardness  traverse-spline  and  body 
o  tensile  properties  -  room  temperature  (RT) 
o  Charpy  impact  transition  data  longitudinal  (LR)  orientation 
o  Charpy  impact  -  (RT  and  -40°C),  transverse  (TR)  orientation 
o  fracture  toughness  (RT  and  -40°C)  (LR  and  TR  orientations) 

Table  2  contains  a  near-complete  chemical  analysis  for  the  GD  material  plus  a 
carbon  and  sulfur  content  for  the  MPU  product.  As  can  be  seen,  the  GD  bar  conforms 
to  the  specifications  for  AISI  8660,  which  is  also  shown  in  Table  2.  It  also  meets 
phosphorus  and  sulfur  requirements  (0.040%  each)  in  MIL-5-45387BAT. 

Table  3  shows  the  results  of  a  Rockwell  C  hardness  (HRC)  traverse  in  both  a 
spline  and  body  location  for  a  GD  material  and  the  spline  location  for  the  MPU  bar. 
It  can  be  stated  that  both  companies  meet  the  hardness  requirements  of  HRC  47-51  per 
MIL-5-45387BAT. 

Tensile  data  from  the  GD  bar  are  shown  in  Table  4.  Duplicate  tests  were  con¬ 
ducted  and  the  data  are  consistent. 

Longitudinal  and  transverse  Charpy  V-notch  energy  data  are  tabulated  in 
Table  5.  Longitudinal  tests  were  conducted  over  the  temperature  range  of  +240°C  to 
-40°C  and  transverse  tests  at  +22°C  and  -40°C.  A  comparison  of  the  data  at  room 
temperature  and  -40°C  shows  that  the  longitudinal  energy  values  are  superior  to 
those  of  the  transverse  orientation.  Values  are  7.7  and  5.5  ft-lb  and  2.7  and 
2.9  ft-lb,  respectively. 

Fracture  toughness  data  (Kq)  are  tabulated  in  Table  6.  These  values  were 
obtained  using  slow  bend  specimens  at  room  temperature  and  -40°C  in  both  an  LR 


and  TR  orientation  from  a  GD  bar  and  room  temperature  TR  orientation  from  an  MPU 
bar.  Room  temperature  LR  values  are  higher  than  TR,  43  versus  38.3  ksi V  in.  GD  TR 
data  is  slightly  higher  than  MPU,  38.3  versus  34.2  ksiVTnT  Two  short  rod  tests 
were  conducted  on  MPU  material  with  resulting  average  Kq  value  of  35.5  ksi V  in 

Microstructural  and  retained  austenite  studies  were  also  conducted  on  the  GD 
material.  The  microstructure  consisted  of  tempered  martensite  and  3%  or  less  of 
retained  austenite. 


IV.  FAILURE  ANALYSIS 

In  addition  to  the  metallurgical  investigations  which  were  performed  on  failed 
torsion  bars,  fractographic  examination  was  also  performed.  One  of  the  basic  prin¬ 
ciples  of  failure  analysis  is  the  preservation  of  the  fracture  surface  in  a  condi¬ 
tion  as  close  as  possible  to  that  at  the  time  of  failure.  Another  principle  is  the 
compilation  of  historic  data  concerning  the  failed  component.  The  first  principle 
issue  was  not  achievable  since  the  available  torsion  bars  were  found  in  scrap  bins 
at  either  the  depot  or  testing  grounds  resulting  in  fracture  surfaces  with  some  rust 
on  them.  The  second  issue  was  resolved  by  using  a  statistical  model  for  the  random 
loading  spectrum,  this  will  be  discussed  later. 

Figures  12a  through  12d  show  typical  field  failures  which  were  obtained  from 
Anniston  Army  Depot,  APG,  Ft.  Knox  and  GD,  which  exhibited  similar  failures.  They 
obviously  initiated  within  the  splines  at  the  ends  of  the  bar,  and  in  most  cases  at 
the  wheel  arm  end.  One  exception  is  a  case  observed  at  APG  where  the  failure  was  at 
the  anchor  block  end.  These  observations  are  supported  by  case  studies  obtained 
from  the  Measurement  and  Analysis  Division  at  APG*  and  GD.^ 

Light  microscopy  and  scanning  electron  microscopy  (SEM)  were  used  in  the  examina¬ 
tion  of  the  fracture  surfaces.  The  failed  bar  in  Figure  12a  was  cut  open  and  examined. 
The  arrow  shows  the  origin  of  the  fracture.  The  arrow  on  Figure  13a  shows  the  surface 
of  the  serration  which  experienced  surface  plastic  flow  at  the  initiation  site. 

This  indicates  high  stresses  at  the  inner  end  of  the  splines.  There  is  a  slight 
indication  of  fatigue  marking  in  Figures  13b  and  13c,  but  not  as  large  as  usual. 

The  SEM  examination  (Figure  13d),  however,  did  not  show  the  usual  striation  beach 
marks  associated  with  fatigue  crack  growth.  This  indicated  that  once  fatigue  crack 
initiation  occurred,  the  propagation  stage  was  very  short.  Further  away  from  the 
initiation  site,  the  fracture  surface  showed  signs  of  corrosion  and  hydrogen 
embrittlement.  It  should  be  noted  here  that  the  fracture  surface  was  old  and  c^uld 
have  been  contaminated  by  the  environment  after  failure. 


V.  STRESS  ANALYSIS 

The  stress  analysis  of  the  splines  was  performed  to  evaluate  various  probable 
modes  of  failure.  For  the  spline  teeth,  the  possible  modes  of  cracking  are  longi¬ 
tudinal,  inclined,  and  transverse,  with  respect  to  the  axis.  Such  cracks  (Figure  14) 
are  associated  with  various  stress  states,  where  the  predominant  ones  are  torsional 

5.  Aberdeen  Proving  Ground  Report,  Case  Nos.:  78-M-4  -  Comparison  Test  of  M60AI  Rise  Tanks:  79-M-9  -  Comparison  Test  of  M60AI 
Rise  Tanks;  80-M-34  •  Failure  Analysis  of  M60A3;  82-M-37  •  Analysis  of  Torsion  Bar  Spring  Failures;  and  82-M-32  -  Failure  Analysis  of 
Two  M60A3  Torsion  Bar  Springs.  Measurement  and  Anal.  Div.,  Phys.  Test  Branch,  STEAP-MT-G,  Aberdeen  Proving  Ground,  MD 

6.  Rebuilt  Tank  Reliability  Analysis.  General  Dynamics,  Land  System  Div.  M60,  MG-8 0-01731-006,  March  1982. 
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shear  stresses,  Hertzian  stresses,  and  bending  stresses.  These  stress  states  are 
three  dimensional  in  nature  and  therefore  the  analysis  could  be  very  expensive  unless 
some  simplifications  are  made.  Three  types  of  stress  models  were  chosen  to  simplify 
the  problem.  The  final  results  were  then  combined  and  used  in  the  failure  evaluation. 

Bearing  Stresses  on  Spline  Teeth 

Due  to  the  existence  of  a  missing  tooth  in  the  spline,  it  is  suspected  that  the 
highest  stressed  area  in  bearing  would  be  in  that  location.  The  bearing  stresses 
were  evaluated  assuming  that  the  load  on  the  missing  tooth  is  distributed  equally  on 
both  neighboring  teeth  (Case  1)  or  only  on  the  forward  tooth  (Case  2).  The  analysis 
was  carried  out  using  sub structuring  and  superposition  of  the  two  load  cases. 

The  stresses  produced  by  the  contact  forces  depend  on  the  number  of  teeth  in 
contact,  and  the  distribution  of  contact  force  along  the  teeth.  The  number  of  teeth 
in  contact  depends  on  the  tolerances,  misalignment  and  load,  while  the  distribution 
of  the  contact  force  depends  on  the  teeth  shape  and  moments  of  inertia  of  the  axle 
and  block.  Figure  15  shows  the  finite  element  model  used.  Figures  16  through  21, 
the  stresses  in  ksi  resulting  from  bearing  forces  are  shown. 

Torsion  of  Nonuniform  Bar 

The  enlargement  of  the  ends  of  the  torsion  bar  leads  to  stress  concentrations 
at  the  fillets.  The  solution  of  this  problem  was  carried  out  using  the  finite  element 
method;  a  heat  transfer  analogy  of  the  classical  Michell  torsion  problem.  Appendix  A 
gives  the  details  of  the  thermal-mechanical  analogy  and  the  boundary  conditions  used 
in  the  analysis.  The  finite  element  mesh  and  the  results  in  ksi  of  this  analysis 
are  shown  in  Figure  22. 

Evaluation  of  Stress  Intensities 

Four  crack  locations  were  considered  in  the  fracture  analysis: 

1.  Surface  crack  in  the  shaft  (Mode  III) 

2.  Subsurface  crack  in  the  shaft  (Mode  III) 

3.  Surface  crack  in  the  spline  area  away  from  the  missing  tooth  (Modes  I  and  III) 

4.  Same  as  3,  but  at  the  missing  tooth  (Mode  I) 

The  Mode  I  analysis  used  the  results  of  the  finite  element  model  in  Figure  15. 

Mode  III  results  in  ksi  of  Cases  1  to  3  were  obtained  by  the  solution  of  a  St.  Venant 
problem  (Figure  23).  The  finite  element  penalty  method  was  used  for  the  multiply 
connected  region  case  of  subsurface  cracks  (see  Appendix  B).  In  all  the  above 
fracture  analyses  the  quarter-point  crack  tip  elements  were  used.^ 

Redesign  -  Reduction  of  Stresses  in  Spline 

As  demonstrated  from  the  failure  analysis  and  from  the  bearing  stress  and  fracture 
analysis,  the  splines  are  the  most  highly  stressed  region  in  the  torsion  bar.  A 
reduction  of  the  contact  stresses  can  be  accomplished  through  a  more  uniform  distri- 


7.  BARSOUM,  R.  S.  On  the  Use  of  isoparametric  Finite  Elements  in  Linear  Fracture  Mechanics.  Int.  J.  Num.  Meth.  Eng.,  v.  10, 
1976,  p.  25-37. 
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bution  of  the  stresses  by  redesigning  the  splines.  Figure  24a  shows  the  bearing 
stress  distribution  from  current  spline  design  (uniform  teeth).  Using  tapered  teeth 
permits  intentional  redistribution  of  the  contact  stresses  between  the  splines  and 
end  blocks.  From  the  deformation  of  the  splines,  as  determined  by  the  finite  element 
analyses  in  section  V  -  Bearing  Stresses  on  Spline  Teeth  and  Torsion  of  Nonuniform 
Bar  -  it  was  calculated  that  a  0.005  in.  tapering  increases  the  contact  zone 
and  reduces  the  maximum  bearing  stresses  by  more  than  25%  as  shown  in  Figure  24b. 
Machining  of  such  tapered  teeth  is  accomplished  using  a  gear  shaper. 


VI.  SIMULATION  PROCESS  AND  RELIABILITY  EVALUATION  DUE  TO  RANDOM  LOADING 


Random  Loading  -  Amplitude  Displacement  Model 


There  is  a  lack  of  load  spectrum  data  for  the  M60  tank.  The  only  available  data 
was  that  of  an  Ml  tank  test  conducted  at  APG  is  shown  in  Figure  25.  The  test  track 
contained  a  series  of  artificial  ramps  and  bumps  (Belgium  blocks),  to  simulate  a 
specific  test  course  (Figure  26). 


The  angular  amplitude  distributions  of  three  bars  from  these  tests  is  shown  in 
Figure  27.  The  positive  and  negative  angular  displacements  of  the  bars  as  a  function 
of  tank  travel  are  shown  in  Figure  27a.  The  zero  angular  rotation  is  referenced  to 
the  static  deflection  of  the  torsion  bar  and  hence  does  not  give  zero  stress  in  the 
bar.  Figure  27b  shows  the  amplitude  distributions  in  a  manner  describing  percent 
time  less  than  a  specific  value  (positive  sign  (+)j  and  percent  time  greater  than  a 
specific  value  [negative  sign  (-)l  ,  e.g.,  25%  leve.  equals  a  -75%  level.  The  positive 
peak  represents  maximum  angular  displacement  uner  load;  the  negative  peak  is  the  max¬ 
imum  rewind  angular  measurement.  The  range  of  angular  rotation  is  defined  as  follows: 


A0  «  6  +  | 6  | 

where  0  *  maximum  negative  angular  displacement,  and 


(1) 


0  -  displacement  from  Figure  27. 


The  Beta  distribution  provided  the  best  representation  of  the  skewed  amplitude 
distribution.  The  dampening  effects  that  occurred  under  load  resulting  from  a  stop 
used  in  preventing  further  angular  twist  of  the  bar  produces  highly  skewed  discrete 
cumulative  probability  values.  The  Beta  function  is  defined  as: 


f(A6) 


r(P+Q) 


r(P)r(q) 

0<A9<1  P,Q>0 


(A0)P_1  (1-A0)Q_1 


(2) 


The  P  and  Q  values  are  selected  in  a  manner  that  provides  the  best  probability 
density  function  (PDF)  for  representing  the  data.  Figure  28  describes  a  typical 
distribution  and  Table  7  shows  the  excellent  fit  between  predicted  (Beta  representa¬ 
tion)  and  actual  test  results.  As  shown  in  the  fatigue  analysis,  section  VI,  Fatigue 
Life  -  S/N  Curve  Analysis,  angles  less  than  20°  represent  stresses  sufficiently  low 
that  infinite  torsion  bar  life  could  be  expected,  therefore,  a  good  representation 
below  this  angle  is  not  essential. 


i 


» . 


Figure  29  shows  a  typical  load  spectrum  result  of  an  M60  A3  torsion  bar.  The 
results  were  obtained  for  a  DADS*  simulation  process  performed  at  TACOM.  A  stochastic 
model  was  introduced  in  order  to  simulate  the  bar  response.  The  spectrum  is  repre¬ 
sented  by  positive  and  negative  angular  displacements  versus  time,  as  shown.  Included 
in  the  figure  are  tabulated  absolute  amplitude  values  using  Equation  1. 

The  uniform  probability  density  function  was  selected  to  represent  the  amplitude 
distributions.  The  selection  process  involved  evaluating  an  informative  quantile  (IQ) 
plot®  of  the  amplitude  data.  In  Figure  30a,  the  IQ  plot  of  data  tabulated  in 
Figure  29  is  shown.  The  straight  line  represents  an  exact  uniform  distribution  with 
the  disjointed  line  representing  actual  data.  Figure  30b  is  a  schematic  of  a  typical 
distribution  model  for  DADS  simulation  data. 

Fatigue  Crack  Growth  Law  for  Estimating  Torsion  b-\r  Life 

Initial  efforts  in  applying  the  Monte  Carlo  method  for  determining  reliability 
versus  cycles  to  failure  of  the  torsion  bar  involved  using  the  crack  propagation  laws. 
The  da/dN  data  requires  extensive  testing,  time,  and  material.  For  this  reason, 
and  due  to  the  short  life  expended  in  crack  propagation  because  of  the  operating 
stress  levels,  it  was  decided  to  rely  on  published  results.  As  shown  in  Figure  31, 
the  percentage  of  life  spent  in  initiating  and  propagating  a  crack  depends  on  the 
stress  level  and  material  toughness.  The  da/dN  relationships  for  materials  metal- 
lurgically  similar  to  the  specified  AISI  8660  material  were  obtained  from  References 
9  through  11,  and  are  shown  in  Figure  32.  The  dry  air  results  made  available  by 
Barsoml®  provided  the  most  representative  estimates  of  crack  growth  versus  stress 
intensity  (AK)  since  the  torsion  bar  is  protected  from  the  environment.  It  is  to  be 
noted  that  the  mean  stress  effect  was  ignored  in  this  analysis.  From  the  basic 
da/dN  relationship,  N  cycles  to  failure  as  a  function  of  crack  growth,  angular  dis¬ 
placement,  and  the  geometry  of  the  region  where  the  crack  initiates  in  the  bar  may  be 
obtained  from  the  following  relationships: 


N 


®f  _ da _ 

ai  0.66  x  10-8  AK2*25 


where  AK=A.A0|/tm 
J 


(3) 


with  A j ' s  (see  below)  representing  the  stress  intensity  multiplier  associated  with 
each  configuration  of  the  cracks  discussed  in  the  previous  sections.  The  A0  is  the 
range  of  angle  of  twist. 

Ai  =  4.91  (missing  spline  tooth) 

A2  =  3.29  (other  spline  regions) 

A3  =  3.26  (torsion  bar  body). 

•DADS  is  a  vehicle  dynamics  computer  simulation  program  used  by  TACOM’s  concepts  lab.  (t  should  be  noted  that  the  authors  are  not 
convinced  that  the  M60  suspension  system  is  accurately  represented  in  the  DADS  program.  Therefore,  the  results  in  Table  9  should 
not  be  taken  as  absolute  values,  but  as  demonstrations  of  the  methodology. 

8.  PARZEN,  E.  Entropy  Interpretation  of  Tests  for  Normality  by  Shapiro-Wilk  Statistics.  Presented  at  the  Twenty-Eighth  Conference  on 
the  Design  of  Experiments  in  Army  Research,  Development  and  Testing,  Monteray,  CA,  20-22  October  1982. 

9.  Damage  Tolerant  Design  Handbook.  Metals  and  Ceramics  Information  Center,  Battelle  Columbus  Laboratories,  MC1C-HB-01,  1975, 
p.  8.2-D. 

10.  Damage  Tolerant  Design  Handbook.  Metals  and  Ceramics  Information  Center,  Battelle  Columbus  Laboratories,  MC1C-HB-01,  1975, 
p.  8.2-F.. 

11.  BARSOM,  J.  M.  Transactions  of  ASME.  J.  of  Eng.  for  Ind.,  Series  B,  v.  93,  no.  4,  November  1971. 
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The  range  m  stress  intensity  factor  AK  was  calculated  by  combining  Modes  I  and 
III  at  locations  1  and  2  on  the  spline  (see  section  V,  Evaluation  of  Stress  Inten¬ 
sities)  through  the  use  of  the  strain  energy  release  rate  definition.  The  fracture 
toughness  in  Mode  III  was  taken  to  be  \/(l  -  v2)/(l  +  v)  multiplied  by  the  Mode  I 
toughness. 

Note,  a  percent  reduction  in  af's  will  provide  a  decrease  in  the  stresses  in  the 
specific  region  of  the  torsion  bar.  The  af  and  the  af  parameters  are  initial  and 
critical  crack  size,  respectively.  The  af  is  obtained  from  the  critical  stress  inten¬ 
sity  value  K^c  for  the  material  considered,  36  ksi  VTn.  in  this  case.  The  angular 
displacement  of  a  bar  A6  can  also  be  represented  by  the  equvalent  stress  range  as 

At  =  rG  (A6)/L  where 
r  *  radius  of  shaft, 

G  **  torsional  modulus,  (4) 

A0  ■  angular  rotation  range,  and 
L  =  length  of  torsion  bar. 

The  Monte  Carlo  Process 

Crack  Propagation  Analysis 

A  schematic  of  the  Monte  Carlo  process  is  outlined  in  Figure  33  for  determina¬ 
tion  of  frequency  of  occurrence  versus  cycles  to  failure  of  the  torsion  bar  using  the 
crack  propagation  law.  An  assumed  normal  distribution  is  used  to  represent  variabil¬ 
ity  in  the  A j ,  af,  and  af  parameters.  A  coefficient  of  variation  (CV)  defined  as: 


establishes  the  standard  deviation  (SD)  for  the  corresponding  known  mean  value  (e.g., 
af  for  initial  crack  size).  Coefficient  of  variation  values  of  5,  10,  and  15  percent 
were  considered  in  developing  the  distributions  in  order  to  examine  the  effects  of 
variability  (inherent  errors  in  measurements,  flaw  size  assumption  or  the  stress 
analysis)  in  the  parameters.  By  selecting  the  above  CV’s,  a  sensitivity  analysis 
can  be  developed,  thereby  providing  a  method  for  recognizing  the  importance  of  the 
parameters  as  related  to  the  cyles  to  failure.  The  A0  Beta  distribution,  as  shown 
in  Figure  33,  has  been  previously  defined  in  Equation  2. 


The  random  numbers  used  in  the  Monte  Carlo  process  are  obtained  from  solving 
for  X  in 


/  fi d* 


where  R  is  a  uniform  random  number  and  ff  corresponds  to  the  desired  type  of  fre¬ 
quency  distribution  for  the  parameter.  A  probability  density  function  for  the  N 
cycles  to  failure  can  be  obtained  by  randomly  selecting  discrete  sets  of  numbers 
from  the  af,  af,  Aj  and  A0  distributions  and  substituting  them  into  Equation  3. 
Note,  there  should  be  an  equal  amount  of  random  numbers  for  each  parameter  to  have 
the  proper  amount  of  numbers  for  the  N  distribution. 


to  v&y:  v to ■>  tovto>  to- 
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Fatigue  Life  -  S/N  Curve  Analysis 


Torsional  bar  life  expectancy  was  obtained  using  the  Monte  Carlo  process  applied 
to  the  S/M  curve  relationship.  The  above  procedure  provided  a  method  for  obtaining 
lifetime  estimates  of  the  bar  by  combining  the  effects  of  crack  initiation  and 
propagation.  A  description  of  the  S/M  curve  is  given  in  Figure  34,  where  the  base 
line  data  was  obtained  from  a  literature  survey  for  material  metallurgically  similar 
to  the  torison  bar  material.  The  survey  provided  a  set  of  S/N  curves  for  torsional 
fatigue,  shown  below,  for  best  representing  the  current  materials  used  in  the  bar. 

Logi0N  *  B  +  0.068  A0 

(7) 

where  B  =  7.70  for  base  line  data. 


The  slope  value  of  0.068  was  essentially  the  same  for  all  curves  in  the  set. 

The  adjustment  in  B  from  7.70  to  8.06  is  made  on  the  basis  of  M60  torsion  bar  quality 
assurance  tests  at  a  single  A0  value  performed  at  the  Scranton  manufacturing  facility 
(see  Figure  6).  A  single  load  equivalent  to  a  42°  angular  displacement  was  applied 
during  the  quality  assurance  torsional  fatigue  test.  Use  of  the  mean  value  and  the 
cycles  to  failure  in  Figure  6  provided  a  more  accurate  estimate  of  B.  The  curves 
representing  a  range  of  10  and  20  percent  reduction  in  bar  stress  are  also  shown 
in  Figure  34. 


The  S/N  curve  Monte  Carlo  process  is  similar  to  the  previously  outlined  method 
for  da/dN  relationship.  The  primary  difference  involves  using  models  for  B  and  A0 
from  Figures  34  and  27,  respectively.  A  schematic  of  the  basic  S/N  representation 
is  shown  in  Figure  35a.  The  simulation  of  S/N  curve  variability  is  shown  for  a 
specific  value  in  Figure  35b.  Figure  35b  also  describes  the  probability  density 
function  (PDF)  for  B.  A  random  selection  of  a  discrete  set  of  numbers  from  A0  and  B 
distributions  is  then  applied  to  Equation  7  in  order  to  obtain  the  Log^oN  value. 

The  process  is  repeated  until  all  values  from  the  two  distributions  are  selected. 
This  process  will  then  provide  a  PDF  to  represent  LogjQN.  Again,  in  the  S/N 
evaluation,  the  mean  stress  effect  was  ignored  for  torsional  loading. 12 


Torsion  Bar  System  Reliability 

By  assuming  a  tank  with  m  torsion  bar  systems,  the  following  procedures  would  be 
applied  in  order  to  establish  the  reliability  of  the  system.  After  the  reliability 
for  each  torsion  bar  failure  is  calculated,  the  next  task  is  to  determine  the  overall 
suspension  system  reliability. 

If  all  of  these  torsion  bar  failures  are  independent  of  each  other  then  the 
system  reliability  Rg  would  be  given  by 


RS  -  n  R;  (8) 

i-1 

where  R£  *  reliability  for  torsion  bar  i,  i  =  1,  2,....,m,  and  m  =  total  number  of 
>  torsion  bars  in  the  suspension  system. 

I 

I 


12.  FORREST.  P.  G.  Fatigue  of  Metals.  Pergamon  Press,  New  York,  1962,  p.  103. 


If  these  torsion  bar  failures  are  not  independent  of  each  other,  then  the  system 
reliability  would  be  given  by 


RS  *  R1  X  8-2/Ri  Rm/Rjjj.jy 


.  • ./Ri 


where 


r"/r»-i/.../r1 

is  the  reliability  of  torsion  bar  m,  given  the  reliabilities  of  m  -  1  ....  1 
bars. 

It  should  be  noted  here  that  the  reliability  of  the  individual  bars  Rj  is 
different  due  to  the  difference  in  load  spectrum  it  experiences. 

Remaining  Life  -  Reliability  of  Operation  After  Specific  Number  of  Cycles 

When  the  torsion  bars  of  a  rebuilt  M60  tank  are  reused  in  service  after  a  Depot 
inspection,  one  would  assume  that  a  specific  number  of  cycles  on  the  bars  could  be 
found.  Currently,  no  information  is  recorded  and,  hence,  the  following  is  a  hypo¬ 
thetical  evaluation.  The  reliability  of  operating  an  additional  number  of  cycles 
when  a  specified  number  of  cycles  of  operation  has  been  completed  is  obtained  in  the 
following  manner.  Initially,  it  is  assumed  that  a  specified  distribution  function 
such  as  f(N)  is  known,  for  example,  the  distribution  of  Log^N  from  the  Monte  Carlo 
method  previously  described.  The  reliability  R(n£,  n)  is  a  conditional  probability 
requiring  the  probability  of  operating  for  (nj  +  n)  cycles  when  n^  cycles  have  been 
completed.  That  is, 


R(nj,n) 


R(n^  +  n) 

R(nj) 


f (N)dN 


n^  +  n 
/“ f(N)dN 


where  n  is  the  additional  mission  in  cycles  after  nj  cycles  of  operation,  the  number 
Ns  (nj,  n)  of  components  (torsion  bars)  that  will  survive  additional  n  cycles  is 
given  by 

Ns  (nj  n)  »  Ng(nj)  •  R(nj  n)  (11) 

where  Ns(nj)  *  number  of  components  starting  the  mission  of  n  additional  cycles. 

Results  of  the  Reliability  Study 

The  proper  number  of  simulations  for  the  Monte  Carlo  method  depended  on  the 
models  under  consideration.  For  example,  5000  and  3000  were  required  for  the  da/dN 
and  S/N  curve  models,  respectively.  Using  a  convergence  rate  criterion  for  the 
calculated  one  percent  values  (Pg  (see  Figure  36)1  and  recognition  of  the  third 
and  fourth  moment  stability  of  the  LogiQN  distribution  provided  an  excellent  method 
for  determining  required  number  of  simulations.  Differences  in  percentile  values 
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for  coefficient  of  variation  (CV)  of  10  and  15  percent  were  minimum  for  da/dN 
simulation  of  parameters.  The  10  percent  value  was  used  for  all  da/dN  calculations. 

The  torsion  bar  reliability  results  from  the  da/dN  relationship  are  given  in 
Figure  37.  The  current  design  results  were  obtained  from  Equation  3,  with  A2  =  3.29, 
which  represents  the  stress  state  in  the  spline  region.  They  indicated  a  relative 
limited  lifetime  range  of  14  to  500  miles,  with  probability  of  survival  (PS)  values  of 
0.99  and  0.01,  respectively.  An  appropriate  increase  in  af  (final  crack  length)  from 
Equation  3  represents  the  40%  increase  in  the  Kjq  value.  This  represents  an  improve¬ 
ment  in  the  material's  capability  with  respect  to  acceptance  of  larger  flaw  sizes 
prior  to  failure.  The  slight  improvement  in  the  bars  capability  indicates  that  an 
improvement  in  material  will  not  significantly  improve  bar  performance.  The  25%  and 
50%  reduction  in  Kj  (stress  intensity)  in  Figure  37  is  obtained  from  reducing  A2  in 
Equation  3  by  the  respective  percentages.  These  reductions  represent  improvements  in 
the  design  of  the  spline  section  of  the  bar  as  shown  in  Figure  24.  The  Kjjj  failure 
in  the  shaft  represent  situations  where  failure  occurs  in  the  shaft  rather  than  in 
the  spline  region. 

The  maximum  life  of  70  miles  at  25  mph  achieved  from  a  50  percent  improvement 
in  spline  design  with  a  0.99  probability  of  survival  indicates  that  there  is  very 
limited  life  of  the  bar  after  crack  initiation.  Table  8  describes  minimum  life 
estimates  (99  percent  survivability)  for  crack,  e.g.,  propagation  of  the  torsion 
bar  with  respect  to  various  tank  velocities  and  the  design  improvements.  Tank 
travel  at  5  mph  (lowest  speed)  with  a  50%  reduction  in  Kj  value  shows  propagation 
life  expectancy  of  only  341  miles  at  0.99  probability  of  survival. 

The  failure  probability  obtained  from  the  S/N  curve  -  Monte  Carlo  application 
is  shown  in  Figure  38.  The  resultant  exponential  form  is  consistent  with  that 
expected  from  the  S/N  model  in  the  analysis. 

A  graphical  display  of  PS  versus  miles  to  failure  is  shown  in  Figure  39  for  the 
25  mph  tank  velocity.  The  life  expectancy  of  the  bar  is  much  greater  than  that 
obtained  from  the  da/dN  analysis.  The  minumum  life  estimate  (0.99  PS)  of  292  miles 
is  21  times  greater  than  the  14  miles  determined  from  the  da/dN  results  for  crack 
propagation.  This  result  indicates  that  most  of  bar  life  occurs  prior  to  crack 
initiation.  Therefore,  the  torsion  bar  should  be  manufactured  in  such  a  manner  that 
flaws  are  minimized.  The  current  shot  peening  used  in  the  manufacture  of  the  bar 
indicated  recognition  of  this  fact  by  the  manufacturer.  The  bar  reliability  estimate 
obtained  after  an  assumed  741  miles  of  tank  travel  are  shown  in  Figure  39.  The  in¬ 
crease  of  the  probability  of  survival  from  0.90  to  0.99  if  the  bar  survives  the 
initial  741  miles  does  not  provide  a  sufficient  gain  to  warrant  reusing  bars  since 
the  minimum  increase  in  expected  life  is  reduced  very  rapidly.  The  results 
(Table  9)  for  a  20  percent  reduction  in  design  stress  gives  865  miles  for  a  PS  of 
0.99.  This  is  a  considerable  improvement  when  compared  to  that  of  292  miles  for  the 
current  design.  The  results  from  velocity  ranging  from  5  mph  to  25  mph  in  increments 
of  5  mph  are  given  in  Table  9  with  respect  to  current  10  and  20  percent  improvements 
in  design.  Reducing  velocities  of  tank  operation  obviously  improves  reliability  of 
the  torsion  bar. 

It  should  be  noted  that  the  field  test  data  used  in  the  comparison  with  the 
reliability  calculation  refers  to  the  failure  of  the  first  (or  front)  torsion  bar. 

The  front  wheel  arms  have  a  different  suspension  system’ than  the  rest  of  the  wheel 
arms,  and  hence  the  spectrum  of  loading  on  the  torsion  bars  is  quite  different.  They 


are  subjected  to  higher  loads  and  have  a  higher  incidence  of  failures.  The  same 
reliability  calculations  can  be  performed  on  the  other  bars  in  order  to  calculate 
the  reliability  of  the  system  (Equations  8  and  9). 

Examination  of  current  design  mileage  capability  of  the  bar  for  20  and  25  mph 
indicates  a  range  from  276  to  292  miles.  These  results  agree  with  the  262-mile 
minimum  life  obtained  from  Aberdeen  Proving  Ground  (APG)  test  results,  Table  10 
(taken  from  Reference  13  of  bar  failure  from  3-mile  test  course),  which  are  repre¬ 
sented  in  the  probability  distribution  of  Figure  40.  This  course  and  tank  velocity 
were  similar  to  those  used  in  obtaining  the  load  spectrum  results.  The  excellent 
agreement  between  the  predicted  and  actual  life  expectancy  of  the  bar  indicates  the 
desirability  of  the  Monte  Carlo  process  for  modelling  variability  of  loads  spectrum 
(design  stress)  and  S/N  curve  (material  capability)  results.  Table  10  gives  the 
failure  results  of  tests  run  at  APG.  The  failure  life  of  the  torsion  bars  are  in 
agreement  with  the  above  predictions. 

Although  excellent  agreement  has  been  obtained,  the  authors  would  have  preferred 
representing  the  spectrum  load  consistent  with  an  individual  peak-to-peak  angular 
displacement.  The  simplification  applied  using  the  negative  peak  as  base  and  repre¬ 
senting  the  displacement  relative  to  this  value  was  a  good  approximation  to  the 
available  individual  displacements.  This  approximation  would  provide  a  slightly 
conservative  estimate  in  the  reliability  values.  Using  the  ASTM  recommended  pratice 
of  representing  the  lower  of  3  times  the  standard  deviation  band  of  the  S/N  curve 
as  a  measure  of  material  fatigue  loading  capability  combined  with  maximum  angular 
displacement  (46  degrees)  for  25  mph,  a  minimum  life  estimate  of  112  miles  resulted 
for  the  bar.  Selecting  this  number  as  a  design  value  could  result  in  an  overly 
conservative  estimate,  since  chance  that  this  maximum  displacement  could  occur  to¬ 
gether  with  the  S/N  curve  as  the  actual  lower  band,  described  above,  is  extremely 
small. 
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A  minimum  life  of  575  miles  was  obtained  by  using  the  maximum  displacement 
value  with  the  original  S/N  curve  where  B  =  8.06.  This  result  is  obviously  wrong 
since  in  the  limited  samples  of  23-bar  failures,  two  of  them  failed  at  mileage  less 
than  400  miles  (see  Figure  40). 

Results  from  the  application  of  the  Monte  Carlo  process  for  DADS  simulation 
data  is  shown  in  Table  11.  The  PS  percentages  and  their  corresponding  miles  of  tank 
travel  are  tabulated  for  road  arms  1L,  2L,  and  6L  at  5,  18,  and  25  miles  per  hour. 
The  APG12  course  (18  mph)  was  the  most  severe  according  to  the  TACOM  representative. 
The  results  from  Table  11  reflect  this,  since  the  0.99  percent  PS  of  14  miles  was 
the  lowest  obtained  from  any  of  the  computed  life  estimates  at  that  PS  value.* 
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Since  this  course  introduces  unusually  large  angular  displacements  of  the  bar 
1L  (see  Figure  29)  and  requirements  are  only  one  failure  in  a  hundred,  this  is  not 
an  unrealistic  estimate  of  bar  life.  It  is  obvious  from  examination  of  Table  11 
that  road  arm  6L  will  survive  much  longer  than  either  1L  or  2L.  Traveling  at  5  mph 
in  such  a  severe  course  will  also  increase  bar  life  considerably. 


! 

A 
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*Sec  footnote  on  page  7. 

13.  Aberdeen  Proving  Ground  Report,  Gtie  No.  82-M-37,  MT-5376,  Aberdeen  Proving  Ground  Test  Branch. 
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VII.  NONDESTRUCTIVE  TESTING  CONSIDERATIONS 

Review  of  torsion  bar  field  failure  reports  and  examination  of  a  number  of 
fractured  bars  clearly  illustrated  that  fatigue  failures  predominately  initiate  in 
the  spline  ends  rather  than  in  the  reduced  cross-section  central  gage  length.  Based 
on  these  data,  it  was  initially  considered  that  a  periodic  magnetic  particle  inspec¬ 
tion  for  fatigue  cracks  in  the  spline  surfaces  would  provide  a  technique  for  pre¬ 
dicting  premature  failures.  However,  results  of  the  subsequent  fracture  mechanics 
analysis  reported  herein  have  identified  the  critical  crack  size  as  0.019  inch. 
Furthermore,  the  analysis  shows  that  an  existing  fatigue  crack,  0.005  inch  deep,  will 
propagate  to  critical  size  within  only  a  few  miles.  Therefore,  if  a  fatigue  crack 
equal  to  or  less  than  0.019  inch  exists,  failure  is  imminent. 

Figures  41a  and  41b  taken  from  Reference  14,  show  that  a  fatigue  crack  in  steel 
must  be  greater  than  0.030  inch  for  a  90Z  probability  of  detection  at  95Z  confidence 
using  the  conventional  magnetic  particle  inspection  method.  Significantly  smaller 
cracks  are  indeed  detectable  using  the  technique,  but  probability  of  detection  will 
be  low  in  a  production  environment,  as  illustrated  in  Figures  41a  and  41b.  Inasmuch 
as  the  torsion  bar  must  be  removed  from  the  vehicle  to  conduct  an  inspection  at  some 
arbitrary  time,  detection  of  a  crack  less  than  0.019  inch  would  be  only  by  chance. 

It  was  therefore  concluded  that,  in  the  case  of  the  M60  torsion  bar,  even  the  most 
sophisticated  nondestructive  crack  detection  scheme  will  not  provide  an  effective 
life  prediction  technique.  Once  a  fatigue  crack  initiates  and  can  be  detected, 
service  life  of  the  torsion  bar  has  essentially  ended. 

VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 

1.  A  methodology  for  obtaining  reliability  of  the  M60  tank  torsion  bar  sub¬ 
jected  to  cyclic  random  loads  has  been  developed  where  probability  of  survival  is 
represented  as  a  function  of  miles  of  tank  travel. 

2.  As  discussed  in  Reference  15,  the  developed  methodology  can  be  applied  to 
other  structures  with  cyclic  random  loads. 

3.  The  use  of  the  method  appears  justified  from  recognition  of  the  excellent 
agreement  between  predicted  reliability  estimates  and  those  obtained  from  the 
actual  bar  life  (miles  to  failure)  experienced  during  the  tank  operation. 

4.  Determination  of  minimum  bar  life  was  21  times  greater  from  the  application 
of  the  S/N  curve  model  than  that  of  the  assumed  da/dN  model.  This  indicates  that 
most  of  the  bar  life  is  expended  in  crack  initiation. 

5.  Application  of  deterministic  procedures  {use  of  lower  three  SD  bound  for  S/N 
curve  (ASTM  method)  and  mean  S/N  curve]  provide  over  and  under  design  value  estimates 
respectively  while  the  Monte  Carlo  method  outlined  in  the  text  accurately  describes 
acceptable  design  values. 


14.  RUMMEL,  W.,  et  al.  Detection  of  Tightly  doted  Flawt  by  Nondestructive  Testing  Methods  In  Steel  and  Titanium.  Martin  Marietta, 
Denver,  CO,  Final  Report  MCR-76-476,  NAS  9-14453,  September  1976. 

15.  BARSOUM,  R.  S.,  et  al.  M60-Tordon  Bar  Reliability.  Army  Symponum  on  Solid  Mechanics,  Newport,  Rl,  October  1984. 
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6.  Calculations  show  that  redesigning  of  the  bars  to  reduce  the  stresses  in¬ 
creases  the  life  considerably.  A  slight  improvement  of  20X  is  much  more  effective 
in  increasing  life  than  requiring  a  higher  toughness  in  the  material. 

7.  Explore  the  use  of  the  torsion  bar  presetting  operation  at  the  depot  as 

a  proof  test  for  reused  bars,  and  follow  it  by  100Z  magnetic  flux  inspection.  Also, 
explore  reshot  peening  of  the  whole  bar,  or  replacement  of  the  pinion  and  anchor 
blocks  at  the  depot. 

8.  Recommend  revisions  to  MIL-S-45387A  by  requiring  the  use  of  aircraft 
quality  material  and  including  requirements  for  grain  size,  cleanliness,  and  Charpy 
V-notch  impact  energy  for  both  room  temperature  and  -40°C.  Such  would 

longer  life. 

9.  For  new  bars,  it  is  recommended  that  materials  other  than  AISI  8660,  such 
as  AISI  4350  steel  be  considered.  Also,  induction  hardening  followed  by  shot 
peening  should  be  explored.  References  16  and  17  cover  this  issue  to  some  extent. 


IX.  TECHNICAL  GAPS 

In  addition  to  the  above  conclusions  and  recommendations  the  study  uncovered 
the  following  “Technical  Gaps"  which  have  to  be  resolved  in  future  studies  of  this 
type. 

1.  Actual  load  spectra  -  the  load  spectra  on  actual  test  courses  and  other 
field  applications  is  not  available.  The  design  of  new  or  improved  suspension 
syste  -ill  require  this  knowledge. 

2.  Data  from  various  manufacturers  -  torsion  fatigue  data  was  only  available 
from  GD.  Other  manufacturers  and  potential  suppliers  of  torsion  bars  had  no  data. 

In  addition,  the  form  of  the  data  generated  needs  to  be  modified  based  on  modern 
statistical  methods  as  discussed  in  his  report. 

3.  In  general,  there  is  a  lack  of  fatigue  and  fatigue  crack  growth  data  for 
materials  used  in  torsion  bars.  In  addition,  environmental,  finishing,  and  mean 
stress  effects  on  fatigue  life  are  very  difficult  to  obtain. 

4.  The  effect  of  shot  peening  and  hardening  on  fracture  and  fatigue  life  needs 
an  extensive  study. 

5.  Currently,  there  are  no  NDE  methods  for  detection  of  small  cracks.  Also, 

NDE  methods  for  field  inspections  are  needed. 

6.  In  order  that  bars  which  exhaust  their  fatigue  life  can  be  retired,  a  method 
of  mileage  record  keeping  for  every  bar  is  needed  (this  could  be  expensive).  Micro¬ 
processors  or  mechanical  angular  displacement  accumulators  would  be  helpful  for  such 
a  task. 


16.  GIBSON,  D.,  BOERMAN,  G.,  and  GORUM,  A.  High  Strength  Tonton  Bar  Development  for  Army  Apptcationt.  Oon tract 
DAAG46-76-C-0074,  FMC  Corp,,  San  Joa e,  CA,  Final  Report,  AMMRC  TR  19-5,  January  1979. 

17.  High  Strength  Steel  Tonkm  Ban,  TACOM,  R7D  Tech.  Report  No.  12746,  30  May  1983. 


7.  Maintenance  procedures  need  to  be  revaluated  in  lieu  of  failure  life  of  the 
component.  In  addition,  the  interaction  between  various  components  should  be  con¬ 
sidered  in  their  retirement. 

8.  Introduction  of  the  censoring  process  should  be  made  in  order  to  reduce  the 
amount  of  fatigue  cycles  necessary  for  QA  of  the  bars.  This  will  reduce  the  cost 
and  time  in  the  testing  process. 
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Figure  1.  M60  tank. 


Figure  2.  Numbers  1  through  6  roadwheel  support  housing,  torsion  bar,  and  related  parts. 


Figure  3.  Torsion  bar  and  teeth  detail. 


Figure  4.  Torsion  bar  endurance  test  fixture  (shown  with  scatter  shield  removed) 


torsion  bar  fatigue  test  results  (1979). 


M60  torsion  bar  fatigue  test  results  (1981), 
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Weibull-censored  data  analysis  of  Scranton  torsional  f 


Weibull-censored  data  analysis  of  Scranton  torsional  fatigue  test  results. 
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Figure  26.  Belgium  block  ramps. 
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spectrum  (run  -  speed  25  mph). 
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Figure  28.  Spectrum  load  representation  (run  no.  48  -  speed  25  mph). 
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Figure  29. 
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Figure  34.  Torsional  fatigue  life. 


Torsional  Fatlaua  Life 


Slaulatioo  of  8  froe  S/M  Curse  Aaprasentatloti 


Cycles  to  Failure  lelatlon 
Log1QK  ■  M  .06*4* 


Probability 

Density 


Cycles  (N)  to  Failure 


S/N  Curve  Log^N  •  8  ♦  068 


(W*  7<*  7*a  »•*  •« 

Figure  36. 


’  -*  .*  V  .**  >  aV 


.*  I  •  .  *  «’  .*  »  ’  .a  '  .  ”  »  \Oe»‘ 

•  .  •  -  •  -  *  »  «  .  *  a  » ^  +  m  m  as  •*  »  W  .  "a  »  a. 


OOCOOOCOCOCOOOCOCOCOOOCOCOCOOO 


Table  1.  SCRANTON  FATIGUE  TEST  RESULTS 


Heat  Code 


Total  Cycles 


Location  of  Break 


Table  1.  Continued 


Heat  Code 

Total  Cycles 

Location  of  Break 

JHL 

140,401 

Test  stopped  per  L. 

R.  suggestion 

JHK 

52,492 

Test  stopped  per  L. 

R.  suggestion 

J8X 

107,417 

62"  from  large  end 

JHP 

222,374 

Test  stopped  per  L. 

R.  suggestion 

JCX 

92,028 

Test  stopped  per  L. 

R.  suggestion 

JHP 

46,249 

Test  stopped  per  L. 

R.  suggestion 

JHP 

135,232 

6“  from  large  end 

JHP 

80,517 

Large  end  radius 

JHP 

214,090 

48"  from  large  end 

JKP 

167,191 

Test  stopped  per  L. 

R.  suggestion 

JKR 

45,092 

Test  stopped  per  L. 

R.  suggestion 

JKS 

298,271 

60"  from  large  end 

JKD 

123,265 

Small  end  radius 

JGX 

75,850 

34"  from  large  end 

JBV 

56,600 

32"  from  large  end 

JU 

136,800 

Test  stopped  per  L. 

R.  suggestion 

JIG 

92,700 

Test  stopped  per  L. 

R.  suggestion 

JLE 

87,250 

Test  stopped  per  L. 

R.  suggestion 

JLF 

231,750 

Small  end  radius 

JLH 

78,440 

Test  stopped  per  L. 

R.  suggestion 

JLW 

52,840 

Test  stopped  per  L. 

R.  suggestion 

JLX 

125,140 

Test  stopped  per  L. 

R.  suggestion 

JLV 

88,500 

Test  stopped  per  L. 

R.  suggestion 

JMD 

180,320 

Test  stopped  per  L. 

R.  suggestion 

3MM 

115,225 

54"  from  large  end 

JNO 

269,811 

Test  stopped  per  L. 

R.  suggestion 

JNA 

217,368 

Machine  failure 

JMP 

267,951 

Test  stopped  per  L. 

R.  suggestion 

JNC 

88,192 

Test  stopped  per  L. 

R.  suggestion 

JNR 

85,715 

Test  stopped  per  L. 

R.  suggestion 

JLT 

95,703 

Test  stopped  per  L. 

R.  suggestion 

JLl 

80,570 

8“  from  large  end 

JU 

175,235 

Test  stopped  per  L. 

R.  suggestion 

JU 

146,881 

Test  stopped  per  L. 

R.  suggestion 

JU 

127,877 

30“  from  large  end 

JU 

244,732 

18"  from  large  end 

JU 

193,882 

Machine  failure 

JU 

227,500 

42"  from  large  end 

JL1 

288,725 

Test  stopped  per  L. 

R.  suggestion 

JOM 

223,629 

Small  end  radius 

JOP 

443,417 

30"  from  small  end 

JKT 

100,607 

60"  from  large  end 

JLT 

49,422 

60"  from  large  end 

JKT 

75,488 

61“  from  large  end 

JKT 

71,488 

Small  end  radius 

JLT 

41,802 

70“  from  large  end 

JLR 

- 

46"  from  large  end 

JOR 

212,084 

Large  end  radius 

JPB 

350,769 

46“  from  large  end 

JOP 

301,591 

24“  from  large  end 

JPB 

303,336 

Small  end  radius 

JPC 

584,366 

18“  from  small  end 

JOS 

268,040 

Machine  failure 

JOS 

228,710 

20"  from  large  end 

JPC 

358,269 

Machine  failure 

JNL 

77,990 

48"  from  large  end 

JNL 

57,409 

44"  from  large  end 

Bar  No. 

New  Cycles  To  Failure 

Location  of 

Break 

A  failed  torsion  bar  from  Ft.  Knox  spline  failure? 
M60-A3.  Mill  be  sent  to  Hickey. 


82 

1 

185,000 

Machine  broke 

82 

2 

148,000 

Machine  broke 

82 

3 

544,000 

Test  stopped 

82 

4 

447,696 

Failed 

82 

5 

240,348 

Failed  (was  broken  in  Ml  machine) 

82 

6 

362,771 

Machine  broke 

82 

7 

367,307 

Test  stopped 

All  New  Spline  Design 

Table  2.  CHEMICAL  COMPOSITION 


Source/Marking 

Weight  Percent 

GD/JN 

MPU  #5399 

Specification  AISI  8660 

Minimum  Maximum 

Carbon 

0.63 

0.54 

0.55 

0.65 

Manganese 

0.87 

- 

0.75 

1.00 

Phosphorus 

0.009 

- 

- 

0.04 

Sulfur 

0.02 

0.02 

- 

0.04 

Silicon 

0.21 

- 

0.20 

0.35 

Nickel 

0.58 

- 

0.04 

0.70 

Chromium 

1.02 

- 

0.04 

0.60 

Molybdenum 

0.15 

- 

0.15 

0.25 

Table  3.  ROCKWELL  C 

HARDNESS 

Source/Marking 

GD/JN 

MPU  #5399 

Location 

Spline  Traverse 

Body  Traverse 

Spline 

HRC 

(48.0  -  50.3) 

(48.5  -  49.9) 

50.0 

Table  4.  TENSILE  PROPERTIES  OF  GD/JN  TORSION  BAR 


Orientation 

0.2%  Y.S. 
(ksi) 

U.T.S. 

(ksi) 

Elon. 

i%) 

R.A. 

it) 

Long. 

231.4 

257.2 

8.3 

28.4 

229.6 

258.1 

8.7 

28.5 

m 3 

2573 

8.5 

253 

Longitudinal  Orientation  (LR) 


Testing  Temp. 

(°C) 

Charpy  Impact  Energy 
(ft-lb) 

Fibrosity 

(%) 

-40  (-40°F) 

5.5 

5 

+22  (+72°F) 

7.7 

10 

+100  (+212°F) 

7.9 

15 

+160  (+320°F) 

8.8 

50 

+180  (+356°F) 

11.0 

80 

+220  (+428°F) 

9.0 

100 

+240  (+464°F) 

8.5 

100 

Transverse  Orientation  ( TR ) 

Testing  Temp. 

Charpy  Impact  Energy 

Fibrosity 

(°C) 

(ft-lb) 

(%) 

+22  (+72°F) 

3.5 

5 

+22  (+72°F) 

1.8 

5 

777 

7 

-40  (-40°F ) 

2.8 

5 

-40  ( -40°F ) 

3.0 

5 

777 

"3 

Table  6.  FRACTURE  TOUGHNESS 

Kq, 

ksi  /in. 

Source/Marking 

Orientation  Testing  Temp.  (°C) 

Table  7.  SPECTRUM  LOAD  (PROFILE  IV  COURSE)  -  BETA  FUNCTION  REPRESENTATION 


Cumulative 

Probability 


Test  Results 
(Degree) 


Beta 

Representation 
0«  (Degree) 


Cumulative  time  probabilities  of  torsional  bar  angular  displacement  0  adjusted  to 
positive  range  by  0A  =  0  +  |0"|  where  6~  =  maximum  negative  angular  displacement. 


Table  8.  MINIMUM  LIFE  ESTIMATES  (99%  SURVIVABILITY) 
da/dN  CURVE  RESULTS 


Velocity 

(mph) 

Mileage  Expected 
(Function  of 

Spl ine  Stress) 
Current  Design 

Mileage  Expected 
(Function  of 

Spline  Stress) 

25%  Reduction 

Mileage  Expected 
(Function  of 
Spline  Stress) 
50%  Reduction 

5 

71.0 

138.0 

341.0 

10 

29.9 

51.9 

143.0 

15 

15.2 

29.3 

72.3 

20 

14.0 

26.9 

66.3 

25 

14.2 

28.4 

70.2 

Table  9.  MONTE  CARLO  RESULTS  FOR  S/N  CURVE  MINIMUM  LITE  ESTIMATE 
_ (99%  PROBABILITY  OF  SURVIVAL)  VERSUS  VELOCITY  (MPH) 


Veioc  rty 
(mph) 


Mileage  Expected 
(Current  Design) 


Mileage  Expected 
(10%  Design 
Improvement ) 


Mileage  Cxpectcd 
(20%  Design 
improvement ) 


69 

74 

9474 

12970 

2C 

00 

3138 

4420 

45 

b33 

1089 

2 

76 

515 

860 

292 

557 

365 

Note:  A  99%  survivability  estimate  of  262  r.nles  was  obtained  from  cumulative  APG 
mileage  on  vehicles  at  time  of  torsion  bar  failure.  Velocity  of  vehicle  during 
tests  was  approximately  15  to  25  mph. 


Table  10.  TORSION  BAR  SPRING  FAILURES  FROM  SELECTED  APG  TESTS,  SEPTEMBER  1979-April  1982* 


Vehicle  TECOM  APG  Report  Number  Vehicle  Position, * 

Type  Project  Number  Number  Failed  Number  and  Side 


Cumulative  APG  Mileage 
on  Vehicle  at  Time  of  Failure 


M60A1E3 

l-VC-080-060-083 

PV-1 

APG-MT-5376 

9 

1R,1L,6L,2R,4R, 

5R,1R,1L,1R 

1643,  1693,  1860,  2387,  2387, 
2387,  2463,  2471,  2686 

li 

PV-2 

II 

3 

1L,1L,6R 

2286,  3086,  3640 

II 

PV-3 

II 

2 

4L,4R 

3025,  3025 

II 

PV-4 

II 

3 

1L.1R.4L 

1127,  2218,  3093 

II 

PV-5 

•1 

5 

1R,6L,5L,1L,6R 

1190,  1988,  3306,  4085,  4171 

ll 

PV-6 

II 

6 

1R,5L,5L,1R,1L,1L 

546,  2217,  2711,  3303,  3834,  3998 

II 

M60Alf 

PV-7 

•1 

4 

1R,1R,1L,6R 

3411,  3641,  3641,  3641 

l-VC-080-060-087 

APG-MT-5299 

2 

1R.6R 

396,  1596 

M60A3 

1-VC-08F-060-033 

APG-MT-5299 

2 

1L,6L 

1402  average 

M60A3** 

l-VC-087 -060-035 

APG-MT- 

2 

1L,1L 

421,  1418 

M48A5**t 

3-WE-100-DIV-008 

APG-MT- 

3 

6L,3R,1R 

1516,  2638,  2659 

*Taken  from  Appendix  II. 12 

^Positions  1,  2,  and  6  have  shocks  attached  to  the  roadwheels. 

^Detailed  failure  report. 5 

‘Discussed  in  greater  detail  in  this  report. 


Table  11.  PROBABILITY  OF  SURVIVAL  VERSUS  MILES  -  TACOM  SIMULATION  COURSE 


Probability 
of  Survival 


5  mph  (33  Displ ) 


Original 


10%  Design 
Improvement 


20%  Design 
Improvement 


14 

16 

1324 

37 

45 

3107 

70 

89 

5032 

76 

238 

9502 

198 

590 

18350 

2424 

3935 

3951 

6245 

7981 

12580 

7246 

11570 

24010 

12490 

20130 

29280 

1523 

2336 

4686 

• 

3618 

5623 

10590 

5717 

8831 

16910 

,*« 

10320 

16020 

29280 

mm 

APPENDIX  A.  TORSIONAL  STRESSES  IN  A  NONUNIFORM  BAR 

Due  to  the  enlargement  of  the  torsion  bars  at  their  ends  there  could  be  high 
stress  concentrations  in  the  filets  caused  by  the  applied  torsional  moments. 

Thermal  analogy  is  used  in  the  analysis  of  torsional  stresses  of  the  nonuniform 
bars.  According  to  the  Michell  torsion  problem, ^  the  angle  of  twist  satisfies  the 
equation 


!_  /rr3  ilv  .  l_  /G  3  Wv 
<Gr  a7}  +  97  (Gr  9z> 


(A-l) 


and  the  shear  stresses  are  given  by 


T  ,  rr  11  T  -  rr  i! 
Tre  Gr  9r*  z0  G  9z 


(A-2) 


where  G  is  the  shear  modulus,  r  is  the  radius  from  the  central  axis  and  Z  is  the 
distance  along  the  axis.  Using  thermal  analogy  of  axisymmetric  heat  condition  and 
setting. 


r,  z  -  z,  T  -  vi,  K 


(A-3) 


where  T  is  the  temperature,  Ky  and  Kz  the  conductivity  coefficients  in  the  Y  and  Z 
directions,  and  q  the  heat  flux,  the  problem  can  be  easily  solved.  By  substituting 
the  above  quantities  in  Equation  A-l,  results  in  the  familiar  heat  conduction 
equation  of  the  axisymmetric  body 


3  /  „  3T.  ,  9  ,  „  3T. 
9y  (yKy  9y)  +  9z  (yKz  9z5 


(A-4) 


It  is  to  be  noted  that  we  could  have  used  the  substitution  Ky  m  Kz  *  Gr3,  resulting 
in  the  plane  heat  conduction  problem.  However,  the  error  in  the  analysis  would  be 
much  greater  due  to  the  combining  of  the  term  Gr3  rather  than  Gr^  in  every  element. 
The  analysis  would  then  be  carried  out  using  the  standard  heat  conduction  finite 
element  codes. 


1.  TIMOSHENKO,  S.  and  GOOD1ER,  J.  N.  Theory  of  Battictty.  McGraw-Hfll,  New  York,  1951. 


APPENDIX  B.  ST.  VENANT  TORSIONAL  STRESSES  OF  MULTIPLY  CONNECTED  REGIONS 


The  existence  of  subsurface  cracks  in  the  torsion  bar  due  to  inclusions  and 
poor  material  is  a  possibility.  The  evaluation  of  stress  intensity  factors  for 
these  cracks  requires  the  solution  of  a  St.  Venant  torsion  problem  of  a  multiply 
connected  region  which  is: 


2  2 

1_£+  .  -2G0 .  (B-l) 

3  3y 

The  stresses  are  determined  from 


T  .  li,  T  ■  -  21 

Txz  3y  yz  3x 

with  G  the  shear  modulus  and  0  the  angle  of  rotation. 


(B-2) 

The  boundary  conditions  are 


o 

3s 

on  the  outer  boundary  and 

jji,  3<|>/3n  ds  -  2G6Ai 


(B-3) 


(B-4) 


where  A^  are  the  areas  of  the  interior  holes  or  crack  regions  and  Tj_  is  the  boundary. 
For  the  case  of  a  shaft  with  one  or  two  subsurface  cracks  or  holes $  can  be  chosen 
to  be  zero  on  the  outer  boundary  and  hence  satisfy  the  boundary  condition 
Equation  B-3.  While  for  other  interior  boundaries  Equation  B-4  can  be  substituted 
equivalently  with  4  *  A£.  The  constants  A£  will  remain  undetermined  and  can  be 
found  from  the  solution  of  Equation  B-l  together  with  the  equilibrium  equations  for 
the  whole  section. 

The  other  alternative  is  to  use  the  condition  3$/3s  ■  0  on  all  boundaries, 
which  also  complicates  the  solution  of  the  differential  Equation  B-l. 

The  approach  used  here  is  based  on  the  penalty  function  method. I  In  the 
equation,  B-l  is  cast  in  a  variational  problem  where 

1  -  /  [<f£>2  +  (f|)2  -2G0*jdv  -  min.  (B-5) 


subject  to  the  condition  3$/3s  -  0  on  all  free  boundaries. 
1.  REDDY,  J.  N.  Penalty  •  Finite  Element  Method t  to  Mechanic e.  ASME,  AMD,  r.  51,  1982. 


Using  the  penalty  function  method,  Equation  B-5  is  written  in  the  following  form: 


ip  -  ' [<i>2  +  >2  -  “<>♦> +  1yp  4  <t)2<is  (B-6) 


where  Yp  is  the  penalty  parameter  and  is  taken  as  a  very  large  number  thus  guaran¬ 
teeing  3<f)/3s  to  be  very  small  and  close  to  zero. 

The  minimization  of  Equation  B-6  gives  the  heat  equations  of  a  two-dimensional 
problem  with  the  same  analogy  as  in  Appendix  A.  The  penalty  members  are  linear  heat 
condition  elements  (beam  like)  with  very  large  (or  infinite)  conductivity. 
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